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Abstract. The nonlinear interaction, due to quantum electrodynamical (QED) 
effects, between photons is investigated using a wave-kinetic description. Starting 
from a coherent wave description, we use the Wigner transform technique to obtain 
a set of wave-kinetic equations, the so called Wigner-Moyal equations. These equa- 
tions are coupled to a background radiation fluid, whose dynamics is determined 
by an acoustic wave equation. In the slowly varying acoustic limit, we analyse the 
resulting system of kinetic equations, and show that they describe instabilities, as 
well as Landau-like damping. The instabilities may lead to break-up and focusing 
of ultra-high intensity multi-beam systems, which in conjunction with the damping 
may result in stationary strong field structures. The results could be of relevance 
for the next generation of laser-plasma systems. 



1. Introduction 

Currently, the development of laser technology and laser-plasma accelerators is 
pushing the limits of the achievable field strengths in laboratories to levels unprece- 
dented in human history [1-5]. The successes in laser-plasma based acceleration may 
even hold the promise of reaching the critical Schwinger limit, when the vacuum 
becomes fully nonlinear [3]. Thus, it is clear that the nonlinear quantum electrody- 
namical (QED) vacuum effects will become important. Another possibility that has 
been pointed out is the formation of plasma channels, evacuated plasma cavities 
which could support immense field strengths. It has been suggested that elastic 
photon-photon scattering could be detected within these system, using the next 
generation of laser-plasma facilities [6,7]. Moreover, a large number of astrophysical 
systems, such as magnetars [8], gives rise to more extreme conditions than one could 
ever produce in earth-based laboratories. As an example of nonlinear QED effects, 
the possibility of photon-photon scattering is perhaps the most prominent [9-11]. 
There has been much interest in this particular effect, both from an experimental 
and an astrophysical point of view (see, e.g., [12-29] and references therein, and [30] 
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for an up-to-date discussion and [31] for an overview). The formulation of this ef- 
fect in terms of the Heisenberg-Euler Lagrangian has been used to investigate such 
diverse topics as photon splitting and magnetic lensing [24-27], direct detection 
via second harmonic generation [18], self-focusing [19], nonlinear wave mixing in 
cavities [21,22] and waveguide propagation [20]. The approach presented here will 
be of importance for both the experimental and theoretical questions that may be 
posed concerning photon-photon scattering. 



2. The governing equations 

Let us here start by giving a short review of the necessary equations. This will serve 
as a guide for the steps to follow. We thus set up a coupled system of nonlinear 
Schrodinger equations (NLSE) for two electromagnetic pulses, and the correspond- 
ing acoustic wave equation for the fluid background. We then apply the Wigner 
transformation to the NLSE, and obtain a coupled set of Vlasov-like equations, 
the Wigner-Moyal equations. In previous work [32,33], it has been found that the 
effect of photon-photon scattering is to introduce modulational and filamcntational 
instabilities in coherent photon systems. Here we will show that the Wigner-Moyal 
system also leads to modulational instabilities, and that these can be constrained 
by Landau-like damping. The results and applications thereof are then discussed. 

The nonlinear self- interaction of photons can be expressed in terms of the Heisenbcrg- 
Euler effective Lagrangian [9-11] 



L = eqF + nel [AF 2 + 7G 2 ] 



(2.1) 

2a 2 h 3 /45mjc 5 w 



where F = (E 2 - c 2 B 2 )/2 and G = cE ■ B. The parameter k 
1.63 x lCT 30 ms 2 /kg represents the inverse of a critical energy density. Here a is 
the fine-structure constant, H is the Planck constant, m e the electron mass, and 
c the velocity of light in vacuum. The Lagrangian (2.1) is valid when there is no 
electron-positron pair creation and the field strength is smaller than the critical 
field, i.e., 

w < m e c 2 /h and \E\ < E„- lt = m e c 2 /e\ c (2.2) 

respectively. Here e is the elementary charge, A c is the Compton wave length, and 
£ crit «10 18 V/m. 

In Ref. [33], the coupled equations 

Acfcoi 
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ck Q2 ■ V | E; 



2fc 2 



V ! E - — - (feoi • V)'j E 1 + — ^-S g + ai S 2 j E 1 = 0, 

(2.3a) 

V 2 - (fc 02 • V) 2 ] E 2 + ^ [^S g + a 2 S^ E 2 = 0, 

(2.3b) 

were derived. They describe the propagation of two electromagnetic pulses E\ and 
E 2 on an incoherent radiation background. Here fc j (j = 1,2) is the unperturbed 
vacuum wave vector (with a hat denoting the corresponding unit vector), A = 14k 
or 8k depending on the photon polarisation state, while 

ai, 2 = 2 - 2fc i • k 02 - (fc i,2 • e 02 a) 2 - [feoi,2 • (k 02 ,i x e 02 ,i)] 2 , (2.4) 
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depends on the relative polarisation and propagation directions of the two pulses 
in vacuum. Moreover, <§ g and §i = eo(\Ei\ 2 } is the energy density of the radiation 
gas and the pulse i, respectively. Here the angular brackets denote the ensemble 
average. 

We note that in the co- linearly propagating stationary case, Eq. (3.5) below yields 
<§ g = 2\§q{§\ + £2), while Eqs. (2.3) may be written as 
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and 



dx 2k i 
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E!=Q, 



E 2 = 0, 



(2.5a) 



(2.5b) 



where e = feoi • fco 2 , and we have chosen the direction of propagation along the 
x-axis. Note that e = ±1 depending on whether the pulses are parallel or anti- 
parallel. For parallel propagating beams, ai_ 2 — 0, and the direct coupling between 
the pulses vanishes. Still, because of the response of the radiation background, the 
pulses are coupled, and Eqs. (2.5) exhibit spatial self-focusing [34]. Thus, we here 
generalise the two-dimensional self-focusing results found perturbatively in Ref. [30] 
and numerically in Ref. [32]. The results of the two-dimensional self-focusing due to 
photon-photon scattering can also be understood in the context of the modulational 
instability exhibited by the NLSE [32]. 

The equations (2.3) are coupled to the acoustic wave equation [33] 



d 2 S g 

dt 2 





+ c 2 V 2 





+ 







2\J §\§2 



(2.6) 



for the radiation gas energy density § g ,\ where (5 
e 2 )(fc 2 • ei). 



ei • e 2 + (fei • fc 2 )ei • e 2 - (fei 



3. Kinetic description 

We assume that the polarisation of the pulses remains constant, and define the 
Wigner functions qj as the Fourier transform of the spatial coherence function of 
Ej,j = 1,2 [35-37] 

Qi(t,r,K) = j^- s J dy cxp(iK-y)(E*(t,r + y/2)E t (t,r-y/2)), (3.1) 

f Note that we have made the split S g —+ So + £g, with S g <g So, in accordance with 
Ref. [33], and transformed away the resulting phase shift term in Eqs. (2.3). 
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where Hk can be viewed as representing the momentum of the individual photons. 
The Wigner function Qj has the property 



\Ej\ 2 ) = J dKQj(t,r,n). 



(3.2) 



The transform (3.1) in conjunction with Eqs. (2.3) leads to the Wigner-Moyal 
equation 

dt 



cfeoj + -, — « - -, — (koj ■ «)feoj 



«0j 



«0j 



Vpj + 2% sin ( | V • V K ) Qj = 0. (3.3) 



where the potentials are defined according to 



Ac/c l ( -S g + -OL\<§2 



and ^2 



Acfc 2 ( -<S g + -oli£\ 



(3.4) 



and the intensities (l-Eil 2 ) are given by Eq. (3.2). Moreover, the sin-operator in Eq. 
(3.3) is defined in terms of its Taylor expansion, the arrows indicate the direction 
of operation, and V K denotes the derivative with respect to k. 

We will now assume that the fields have perpendicular polarisations, i.e. /? = 0, 
and Q.j = or 4 if the beams are co- or counter-propagating, respectively. Equation 
(2.6) then simplifies considerably, and takes the form 



dt 2 



— ^-V 2 «?„ = — — XSq 



dt 2 



+ c 2 V 2 (A+^2), 



(3.5) 



4. Stability analysis 

In order to analyse the stability of the system (3.3) and (3.5), we linearise according 
the following scheme. Let Qj = Qoj(K) + Qj(n) exp(iK ■ r — iQt) (and the corre- 
sponding expression for <£*,), where Qj <C Qoj, and <§ g — § g exp{iK ■ r — iQt). For 
co-propagating waves, Eqs. (3.3) (in the Vlasov limit) and (3.5) give the dispersion 
relation 

1 = Wilt + W 2 4, (4.1) 
while in the case of counter-propagating pulses, we obtain 

(§A^ A) 2 (1 - - W2I2) = (1 + |A^ A) WiWiI+l; (4.2) 

where Wj = (4/9)A 2 c/c 0j A, A = (Q 2 + c 2 K 2 )/{Q 2 - c 2 K 2 /3), and 

j± = f j k K ■ V ^gQj u o-) 

i - J T cK x (l- K x /k 0j ) - (c/k 0j )K • k ' { - 0) 

where K x = k n ■ K For q 2 = 0, we obtain a result similar to that of Ref. [44] 
concerning Landau damping. 

For mono-energetic beams, i.e., Qoj(^) = (\Eoj \ 2 }S(n — &oj), the integral (4.3) 
can be reduced, and we obtain the well-known beam modulational instability. In 
the case of parallel propagating beams, with Eqi = -E02 = Eq, koi = fco2 = ko, and 
«oi = K 02 = k n x, we obtain the growth rate [33] 
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where the dimensionless parameters are defined according to 770 = A<£b and rj p = 
XS'p. Furthermore, J? = cK x + iF and K± = K — K x x. 

In general, however, the background distribution is not mono-energetic, and due 
to the spectral width of go we will have Landau- like damping manifested by the poles 
of the integral (4.3). In the one-dimensional case, we can take the spectral width 
kw into account by considering two equal incoherent backgrounds with qo{k) — 
{\E \ 2 } I (V2ttk w ) exp[- k, 2 /2k 2 v ], so that 1+ = 0, and 



r = 



\E \ 2 )k 



1 



iy/nk (O + cK) 



+1 



2ck w {' 2\f2cn w K 
^{\E \ 2 )k 2 (n + cK) 
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exp 
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8c 2 k 2 w K 2 



Ac 2 n z w 



V2K 



exp 



k 2 (n + cK) 



2 1 



8c 2 n 2 w K 2 



(4.5) 



Here we clearly see that the effect of the non-zero spectral width is to introduce a 
damping. 



5. Discussion and conclusion 

The existence of a modulational instability for mono-energetic beams strongly sug- 
gests that the full nonlinear effect due to photon-photon scattering of incoherent 
waves in a radiation background should be taken into account, since any small per- 
turbation can grow to form a large amplitude structure, given sufficient time. On the 
other hand, as the example with a one-dimensional Gaussian distribution shows, we 
must also expect the evolution to be damped by the resonant interaction between 
incoherent modes and fluid modes. Thus, it is not unlikely that the full system can 
account for very interesting structure formation, where the initial growth is gov- 
erned by a modulational type of instability, and that this growth is stabilised by the 
Landau-like damping at a later stage of the dynamical phase of the photon system. 
Thus, one can therefore conjecture about the existence of three-dimensional stable 
photon structures generated by the vacuum nonlincaritics. In order to investigate 
this conjecture, a numerical analysis of the full system would be necessary. 

Situations where the effects presented in this paper may occur range from earth- 
based laboratory systems, such as ultra- high intensity lasers [1,2] and plasma ac- 
celerators [3,7], to astrophysical scenarios, such as the early Universe [30], gamma 
ray bursts [45] and magnetars [8]. If the presents effects do occur, they could also 
show up on small angular scales within high precision cosmology measurements, 
such as the ones presented by WMAP [46-48] . A common feature for all these sys- 
tems is that they at some stage of their evolution lead to extreme radiation energy 
densities, which is a key feature for probing photon-photon scattering. Thus, there 
are good possibilities that in the near future one can put QED through new tests, 
both in laboratory and astrophysical environments, where the latter is perhaps the 
most fascinating, since it will connect very small scales to very large scales as a 
fundamental theory test-bed. 

In the present paper we have used a system of nonlinear Wigner-Moyal equa- 
tions. These equations are coupled through a radiation fluid background, where the 
dynamics is determined by an acoustic wave equation driven by the incoherent pho- 
tons. Using these equations, we have shown that incoherent electromagnetic pulses 
can transfer energy between each other by means of the radiation fluid background. 
Consistent with previous work, we have moreover shown that the resulting system of 
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equations is subject to a modulational instability, and that a nonzero spectral width 
of the background incoherent photons gives rise to a damping, much like the well- 
known Landau-damping of electrostatic waves in plasmas. The implications of the 
results have been discussed, and it was conjectured that stable three-dimensional 
electromagnetic structures may form as a result of photon-photon scattering. 
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